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parameter, , in terms of the stability parameter z/L or the Richard- 
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I. INTRODUCTION 



The advent of more complex and complete numerical models for analyz- 
ing and predicting large scale atmospheric motion has also led to better 
specifications of the boundary layer. The boundary layer is that region 
in the first kilometer over the sea which is defined principally by turbu- 
lent transfer of momentum and heat. Progress has been made in specifica- 
tions of mean surface wind fields and hydrostatic stability conditions. 
These are important parameters for predicting changes in the upper part 
of the ocean or for estimating boundary fluxes in medium range atmospheric 
prediction models. 

The Navy, however, also requires surface layer specifications with 
regard to their effect on the propagation of optical waves, and the 
Naval Postgraduate School has been tasked with providing the basis for 
prediction of the optical properties of the atmosphere for relevant atmos- 
pheric window wavelengths. The objective of such a program would be to 
relate the optical properties and the optically relevant meteorological 
parameters of the marine boundary layer to the observable bulk atmospheric 
parameters . 

The micrometeoro logical quantities which define the bulk meteorologi- 
cal quantities are time-averaged values of wind, temperature and relative 
humidity measured at multiple levels within the constant flux layer. 
Optically relevant parameters measured are turbulence induced temperature 
and humidity fluctuations, their frequency spectra, and the aerosol and 
fog concentrations. The optical properties measured are the modulation 
transfer function of the atmosphere, the image wander, its frequency 
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spectrum, the extinction due to absorption and aerosol/fog scattering, 

2 

and the optical turbulence-structure function, C , related to scintilla- 

n 

tion and defined in more detail in section II 0 It becomes of prime im- 
portance in the investigation of optical properties in a marine 

2 

environment to compute easily • Fortunately, a relationship exists 

2 2 
between and the temperature- structure function, , which is derived 

from micrometeorological considerations. 

A serious handicap in describing the marine environment in the past 
has been the lack of data and the fact that available expressions, which 
are empirical in nature, have only been validated over land. 

This study is an examination of data collected over a period of two 
years under open ocean conditions from the R/V ACANIA. The data are com- 
puted with expressions derived by Wyngaard et al (1971) and Friehe (1976) . 

2 

These expressions relate C T to a stability parameter, either the 
Richardson number or z/L. Both the temperature-structure function and 
the stability parameters are defined later in this manuscript. 
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II • THEORETICAL BACKGROUND 



A. GENERAL 

Turbulence properties of interest in optical propagation are those 
which could be used to describe the intensity of fluctuations in the 
refractive index and that which could be used to describe the size or 
scale of the refractive index inhomogeneities. The former has been re- 
lated empirically to scintillation, beam spread and beam wander; and 
the latter along with the former to the image resolution which is im- 
portant with respect to seeing conditions. 

On the basis of the isotropic nature of small scale fluctuation, 
only one parameter is necessary to describe the intensity of the refrac- 
tive index fluctuation over many scales. It is the refractive index- 
structure function parameter, previously mentioned, 

2 2 2/3 

C = [n (x) - n (x+r ) ] /r 7 (1) 

n 

where n(x) and n(x+r) are refractive indices at two points on a line 

oriented normal to the mean wind direction and separated by the distance 

r which is less than the outer scale, L , (lower end of the inertial 

o 

subrange) and greater than the inner scale, SL (associated with the 

upper limit of the inertial subrange) . The inertial subrange in this 

spectra is the wave number region where energy is transferred without 

dissipation from lower to higher wave numbers. See Figure 1. 

The refractive index is determined primarily by density fluctuations 

2 

and can, therefore, be related to temperature fluctuations. is 
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Gain 




Figure 1. A schematic drawing showing spectral transfer leading to 
existence of the inertial subrange (Lumley and Panofsky, 
1964 ) . 
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2 

related to the temperature-structure function parameter, C , by the 



relation: 



2 .79 x 10 



-6 




C 



n 



(2) 



where C 



2 



can be defined by the temperature variance spectrum, below. 



T 



This turbulence parameter is the most important variable in determining 
the optical propagation characteristics of the atmosphere and will be 
discussed in the next section. 

The following expression for the one-dimensional spectral density 
representation of temperature fluctuation is derived from similarity 
theory and dimensional analysis: 



As described earlier, turbulence is nearly synonymous with tempera- 
ture fluctuations; hence it is necessary to describe the mean thermal 
stratification. Parameters for this are the Richardson number defined 
by profile measurements or (z/L) defined by flux measurements. These 
concepts will also be expanded in a later section. 

B. CALCULATION OF THE TEMPERATURE- STRUCTURE PARAMETER 

2 

There are several methods available to estimate C^ directly and in 

directly. These have been partially introduced in previous discussions 

2 

The first approach is used to determine C ^ directly by the 2/3-law 
expression for the structure function. In a locally isotropic field, 
the temperature-structure function has the form 



S (k) = 0.25 C 2 k 5/3 
T T 



(3) 



C, 



2 



[T (x) - T (x+r) ] /r 



2 2/3 



(4) 



T 
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at separation r of the order of inertial subrange scales. 

A second approach is based on a one-dimensional temperature spec- 
trum, S^Ck), which by definition is the Fourier Transform! of the corre- 
lation function with separation r in the streamwise direction: 

C T 2 = 4 S T (k) k 5/3 (5) 

where k is the streamwise component of wavenumber. Since velocity and 

temperature fluctuations are measured at a fixed point in the flow, the 

resultant spectra are realized at a temporal frequency, f. To obtain 
2 

C T , the temporal (f) and space (k) scales are assumed to be related by 
Taylor*s frozen-field hypothesis, i.e., 

k = 2uf/U (6) 



where U is the mean wind speed at the measurement level. "Frozen” 
implies that the turbulence pattern remains unchanged as it sweeps past 
the probe. 

The third approach involves measuring separately the rate of dissipa- 
tion of turbulent kinetic energy (£) and the rate of dissipation of 
temperature variance (N) and applying the results to an expression 
suggested by Corrsin (1951) , as another form of the temperature spectrum, 
viz. , 



s T (k) = Bne 1/3 k" 5/3 



(7a) 



which yields the following expression for based on its definition in 

Equation (5) : 



C T 2 - bn,- 1 / 3 



(7b) 
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(3 is an empirical constant evaluated to be 0.25). From this last 

2 

form, indirect estimates of C can be obtained based on mean condi- 

T 

tions. This is because £ and H are easily related to boundary fluxes 

and profiles, if steady and horizontally homogeneous conditions exist. 

2 

Expressions which relate to mean properties of the boundary layer, 

such as Ri or z/L, are desirable because the small scale measurements 

are impractical to obtain in most operational or tactical situations. 

2 

Expressions relating to £ and N, Equation (7a) , and similarity 

theory predictions for the dependence of £ and N on momentum and heat 
fluxes were obtained from examinations by Wyngaard et al (1971) : 



u * z 



= f x (z/L) 



N 

T * u * z 



= f 2 (z/L) 



( 8 ) 

(9) 



1/2 3 

where U. = (-u'w 1 ) , T. = -w'T'/U., L = -T U. /kg w'T' and f , (z/L) 

X X a Q rC 

and f^(z/L) are empirically determined functions. 

Direct substitution of Equations (8) and (9) into Equation (7a) 
yields an expression of the form 

C T 2 = T* 2 z" 2/3 f 3 (z/L) (10) 

where f^(z/L) results from the combination of f^(z/L) and f^(z/L). 
Furthermore, since z/L and Ri can be functionally related, a parallel 
dependence on Ri can be obtained, which is (Wyngaard et al) 

C T 2 = z 4/3 0t/ 32) 2 f 4 (Ri) (11) 
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This final expression provides the desired dependence of on 

more readily measured parameters (z, 3T/3Z, and Ri) . 

C. CONSIDERATION OF THE RICHARDSON NUMBER 
1. Derivation of the Richardson Number 

McIntosh and Thom (1972) provide a clear and concise explana- 
tion of the Richardson number which is the ratio of the buoyant kinetic 

energy (K.E.), to the natural kinetic energy of turbulence (K.E . ) of 
b n 

a fluid particle in a turbulent boundary layer. 

Two reference levels (z^ and z^) are chosen within a turbulent 
boundary layer (shown in typical profiles of velocity, temperature and 
density in Figure 2) , such that, z^ = + Z, where Z is the local mix- 

ing length. In other words, turbulent exchange of each property occurs 
through the same characteristic distance, or mixing length, Z. 

The total kinetic energy of a fluid particle of volume V within 
this boundary layer is made up of two independent components: one asso- 
ciated with the horizontal progress of the particle as part of the mean 
flow, and the other associated with its own random movements as part of 
the forced turbulence present in the boundary layer. (The latter is 
related to the work required to raise the particle from z to z^ by 
buoyant action and to the kinetic energy associated with the purely 
vertical movements of the particle.) 

The instantaneous velocity of the particle exceeds the mean value 
by an amount U 1 given by U(z 2 ) - Ufz^); i.e., to a first approximation, 

U f = £(3u/3Z). McIntosh and Thom also assume that it is justifiable 
to equate U* and w 1 , so that the vertical component of the 'natural 1 



18 




19 



Figure 2. Typical profiles of velocity, temperature and density (McIntosh and Thom, 1969) * 



kinetic energy of turbulence, (K.E . ) ^ can be written 

(K.E.) = ^ pVw' 2 = ^ pVU' 2 = i pVil 2 (aU/3z) 2 (12) 

n z z z 

Due to turbulent exchange, the fluid particle initially at z^ 
and with density appropriate to that level, will find itself at 
level z 2 and in surroundings of density p^. At level z^, therefore, 
the particle is subject to a positive buoyancy force of magnitude 
g(p 2 - p^)V = Assuming that a linear temperature profile exists 

locally between z and z 2 and that the buoyancy force acting on the 
particle increases uniformly from zero at z^ to F at z 2 , the work done 
(W) by buoyant action is equal to the product of the force ~ F and the 
distance (z 2 - z^) : 

w = \ g(p 2 - P x ) (z 2 - z x )V (13) 



This amount of work goes directly to increase the upward 
kinetic energy of the fluid particle, equal to (K.E.)^, by an amount 
(K.E.)^' to which W may be equated: 



(R.E . ) b = 2 g (p i - P 2 )(z 2 - z l)V 



(14) 



Finally, on substitution of £ for (z 2 - z^) , £ -g— for (p 2 - p^) 

and - ~ for -j—-, assuming an incompressible fluid, the equation becomes 
T dZ dZ 



(K. Eo ) b = -i P V£ 2 f|| 



(15) 



The magnitude of the ratio (K.E. ) . / (K.E. ) indicates the rela- 

b n 

tive importance of free and forced convection in determining the 
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structure of a fluid boundary layer (shape of velocity profile, nature 
of turbulence, etc.). The negative of this ratio, or the Richardson 
number Ri, is a measure of the stability or instability in a region of 
boundary-layer flow. In lapse conditions Ri is negative, and in inver- 
sion conditions, it is positive. After cancelling appropriate terms 



g 8T/3Z 

T (3u/3z) 2 



(16) 



which is one of the forms used in this study 0 

2 . Other Relationships for Richardson Number 

The Richardson number Ri can also be expressed in terms of 
(z/L) or as a bulk aerodynamic formula. 

The stability of the surface layer regime determines flux in- 
tensities. Monin-Obukhov theory determines the scaling level, L, which 
is proportional to the level where mechanical and thermal production of 
turbulent kinetic energy are equal. The ratio of the height of the 
measurement (z) to the Monin-Obukhov length (L) serves as a stability 
index (z/L) . The Monin-Obukhov length is defined as 



T V U * 

q k t* 

v 



(17) 



where T* = T* + 0.61 q* T, q* = -w'q'/U^, and K = 0.35 (von Karrnan 
v 

constant) . 

In unstable conditions the stability parameter z/L is approxi- 
mately equal to the Richardson number Ri . The following relationships 
between the stability parameter and Ri have been proposed by Dyer and 
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Hicks (1970) for an unstable 


case and Webb 


(1970) for the stable 


case 


z/L = 


Ri 


(stable) 


(18) 


1 - a Ri 


z/L = 


Ri 


(unstable) 


(19) 



Here a is an empirically derived constant equal to 0.5. 

A bulk aerodynamic formula is one which relates fluxes across 
the air/sea boundary to the wind speed at a level and the temperature 
difference between that level and the surface. The derivation of such 
a formula involves several assumptions regarding the stability condi- 
tions of the boundary layer and the turbulent processes within it. If 
valid, however, this type of formula is very useful for most practical 
needs . 

Friehe (1976) selected bulk aerodynamic formulae for empirical 

parameterizations of the surface fluxes of momentum, sensible heat and 

moisture in terms of average surface layer variables. He used the 

selected bulk aerodynamic formulae to derive an expression for the 
2 . 

estimation of in terms of the relatively easily-measured quantities 

mean wind, U, temperature, T, and water vapor density, 0. This formu- 
lation by Friehe is presented in the following section. 

D. EXISTING EXPRESSIONS AND RESULTS 

Wyngaard et al (1971) examined forms of nondimensional expressions 

2 

relating atmospheric stability to the structure parameter, C T . Dimen- 
sional analysis of the height (z) , the friction velocity (U*) , the 

- 1 2 

mean potential temperature gradient (30/9Z) , e, N and produce 

"^In this study virtual potential temperature gradients were used 
to account for humidity gradients over open ocean conditions. 
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functional expressions dependent only on the Richardson number. Simi- 



larity principals are elucidated by using the following Richardson 
number definition to indirectly obtain a comparable temperature-struc- 
ture parameter, i.e.. 



Ri 



g/T Q0/3Z) 

(3u/3z) 2 



( 20 ) 



This relationship differs from Equation (16) in that the temperature 
parameter here is a mean potential temperature. When expressed in 
the finite difference form, a simple approximation to 3U/3Z, for 
example, is 



3U y AU 
3z z^ Az 



U(z 2 ) - U(z ) 
Z 3 S ’ n(z 2 /Z l ) 



(21) 



where the computed number applies to a geometric mean height defined 

. . . 1/2 

by z 3 = ( Zl z 2 ) . 

Equations including the Richardson number functions (f ^ , f^ and 
f^) are equations (8) , (9) , and (11) . Equation (11) relates the tem- 

perature-structure parameter and stability parameter and allows for an 
empirical determination of the £ ^ function at varying conditions of 
stability. If the function f is defined, the temperature-structure 
parameter can be readily calculated from more easily measured para- 
meters, 30/3Z and Ri. In Figure 3 the predicted curve and the data 
obtained from the extensive 1968 AFCRL study of turbulence structure 
over a flat, unobstructed Kansas plain are presented. As will be shown 
in the section describing results, available marine data does not appear 
to agree with this predicted curve. 
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(o©/az) 




Ri 



Figure 3. The dimensionless temperature-structure parameter vs. the 
Richardson number (Wyngaard et al, 1971) . 



2 *. 



Friehe (1976) combined the semi -empirical equation of Wyngaard et 

al with the bulk aerodynamic formulae to provide alternate equations 

2 

for the estimation of in the atmospheric boundary layer above the 

2 

sea. These equations relate to measurable quantities of average 

wind speed and air/sea temperature and moisture differences. 

The expression derived by Wyngaard et al and considered by Friehe 
is Equation (10) with f^ defined as 

-2/3 

f^ = 4.9(1 - 7(z/L)) for z/L <_ 0 

f = 4.9(1 + 2.75 (z/L) ) for z/L j> 0 

and z = height, m 

T* = characteristic temperature scale, K 
The bulk aerodynamic formulae considered by Friehe were 



where 



u * 

w0 

wq 

U 



c 

c 

c 



D 

H 

E 

h 



C D 


5 2 






(22) 


C„ 


U (T 


- T ) 


= C UAT 


(23) 


H 


s 


a 


H 




C E 


0(Q s 




= c uAo 

E 


(24) 



average wind speed at height h, m/s 

average potential air temperature at height h, K 

3 

average water vapor density at height h, g/m 
average sea surface temperature 

average water vapor density adjacent to the sea surface, 
usually obtained by assuming the air to be saturated 
i.e., at dew point T , g/m^ 

drag coefficient, from experiment 

sensible heat flux coefficient, from experiment 

moisture flux coefficient, from experiment 

reference height, usually 10m above the average sea 
surface 
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w = vertical velocity fluctuation, m/s 

3 

q = moisture fluctuation, g/m 
0 = potential temperature fluctuation, K 



The resulting expressions derived by Friehe and adapted for this 

study were the following: 

For -20 < UAT < 25 (mK/s) : 

-3 -2/3 

3.12 x 10 (1 + 1.62X) 7 unstable 

C 2 z 2 ' 3 

_T (25) 

(AT) 2 3.12 x 10 -3 (1 - 0.635X) stable 

•7 At 2 

where X = =— (1 + 0.212 R AQ/AT) Ks /m 

U 

For UAT >25 (mK/s) : 

C 2 z 2/3 

^T Z -3 -2/3 

— r = 8.03 x 10 (1 + 2.60X’) ' (26) 

(At) 

where X’ = (1 + 0.132 R 4 AQ/AT) 

U 



and are the ratios of the coefficients used in making adjust- 
ments to the levels used for this presentation. Results obtained by 
Friehe (1976) appear in Figure 4. 
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Figure 4. The dimensionless temperature-structure parameter vs. X (Friehe, 1976). 



III. THE EXPERIMENT 



A. THE PLATFORM AND LOCATION 

Observations were obtained aboard the R/V ACANIA anchored off 
Monterey, California in Monterey Bay. Measurements were made at four 
levels on two masts spatially separated on the forward deck of the 
ship. An illustration of the sensor locations of these masts is shown 
in Figure 5. The vane and probe arrangement appears in Figure 6. The 
vane maintained a one-dimensional profile of the wind during measure- 
ments . 

The location of Monterey Bay provided an ideal site for the experi- 
mental program as shown in Figure 7. Open ocean differs from land in 
the effects of wave action on turbulence, in the nature of the aerosols 
and fog, and to a lesser degree but perhaps significantly on the humid- 
ity fluctuations. These conditions can, of course, be best obtained 
far at sea. The cost of such activities makes it desirable to work 
near land. Pt. Pinos and Monterey Bay provide very nearly the ideal 
situation. 

Pt. Pinos projects northwest from the mainland toward the prevailing 
northwest wind. Even under storm wind conditions, the wind comes from 
the southwest, still bringing sea air toward shore. The ocean depth sur- 
rounding Pt. Pinos is shallow enough that the ACANIA can anchor at any 
range up to about 15 km to the north and northeast. 

The effects of land in producing acceleration of the airflow are 
very small in the Pt. Pinos-Monterey environment. The shoreline south- 
east from Pt. Pinos is parallel to the prevailing wind. It rises 
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legend 



U Cup Anemometer 
T Quart/; Thermometer 
^ Humidio meter 
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Figure 5. Mounting arrangements aboard the ACANIA. 
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Figure 6. Vane and Probe Arrangement 
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slowly along the water's edge. The eastern edge of Monterey Bay con- 
sists of stabilized sand dunes. The sand beach shore is the mouth of 
the Salinas Valley, about 20 miles wide at the shore line. This valley 
extends approximately 100 miles to the southeast. In daylight hours, 
it adds a convective airflow to the southeast. The conventional day- 
night sea breeze reversal often roughly cancels the prevailing airflow 
at night to produce a calm period rather than a reversal, so that a 
relatively small amount of land air is transported out over the water. 

B . INSTRUMENTATION 

1 . Mean Measurement Systems 

The mean wind measurements were made with a Thornthwaite Asso- 
ciates cup anemometer wind profile register system, model number 104. 

In operation the shaft of a three cup anemometer unit serves as the 
shutter between a light source and a photocell for each revolution. 

The cups are plastic cones reinforced with aluminum frames. They are 
attached to the rotating shaft by stainless steel tubing spaced at 
120 degree intervals about the shaft as shown in Figure 8. The three 
cup assembly sets along with the other sensors were positioned at four 
levels on the bow with electrical leads to the after deck house labora- 
tory. The sets have the characteristics of low starting speeds with 
a small amount of internal friction which aids in checking inertial 
overshooto 

Temperature sensitive quartz crystal probes, (Hewlett Packard 
model HP-2850) were used to measure mean temperatures at each level. 

RF signals from the crystal probes and from a reference oscillator were 
mixed in the HP-2801A readout unit to produce a beat frequency whose 
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Figure 8. C. W. Thornthwaite Anemometer Cups. 
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signature can be analyzed to within 0.001 degrees centigrade per hertz. 

Each sensor simultaneously received pre-experiment calibration against 
a platinum resistance wire thermometer in a temperature controlled cir- 
culating water bath over the expected temperature range. The accuracy 
in achieving a 0.005 degree centrigrade correction factor was a constant 
for each probe. A 3.7 meter flexible coaxial cable is permanently 
attached to the sensor head and the mast mounted probes are housed in 
an aspirated shelter as depicted in Figure 9. Temperature values were 
automatically recorded on a printer tape. 

The Hygrodynamics Digital I hygrometer indicator and the Dunmore- 
lithium chloride sensor were used to measure the relative humidity. The 
Dunmore-type lithium chloride sensor is also pictured in Figure 9. The 
equipment operates on the basis of resistance change within an electro- 
lytic solution generating a reference voltage variance which is propor- 
tional to the relative humidity change. Automatic temperature compensation 
in the instruments meets the following specifications for relative humidity: 

+_ 3% relative humidity below 90% 

Hh 4% relative humidity above 90% 

Calibration of the sensors was accomplished by comparative methods using 
saturated salt solutions in an enclosed chamber. Relative humidity was 
observed and recorded similar to that of temperature as a printed output. 

2. Temperature Fluctuation Sensor 

The temperature fluctuations were measured using a bridge 
developed by personnel at GTE Sylvania, the GTE Sylvania Model 140. 

The system was slightly modified for use in this study. 

This system is a small temperature difference sensing device 
originally designed to be carried aloft by a balloon and equipped to 
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Dunmore-type Lithium 
Chloride Sensor 



C.C. Breidert Company 
Air-X-Houster Type 6L 






Quartz Thermometer Probe 



Figure 9. Photographs of aspirated shelter, hygrometer and 
thermometer probe. 
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process and transmit, by pulse-rate modulation on a radio-frequency 
carrier, this information to a ground-based receiving station for 
demodulation and recording* 

The baseband portion of this system is basically a balanced 
Wheatstone bridge excited by a 3 kHz signal with a synchronous detector 
on the output. Segments of a very small diameter platinum wire serve 
as temperature sensors in opposite arms of the bridge. In the single 
wire mode, one wire is replaced by fixed resistance. The resistance- 
temperature coefficients result in an output from the bridge which is 
proportional to the temperature difference between the two probes or 
a temperature change of a single wire in the single wire mode. 

The sensor wire is 0.5 cm long and 2.5 x 10 ^ m in diameter. 

This extremely small mass allows a response to temperature variations 
of up to 1 kHz while electronic amplification allows temperature differ- 
ences as small as 0.004 degrees C to be observed. 

C. ANALYSIS PROCEDURES 

This particular study was divided into two phases: 1) individual 

mean wind U, mean temperature T, and mean humidity q values were plotted 
on logarithmic scale. Best fit lines were drawn to the points from 
which data were picked off and applied to the expressions developed by 
Wyngaard et al; and 2) values of U, T, and q from individual levels were 
applied to the bulk aerodynamic expressions examined by Friehe. 

Fluctuating temperature T* and profiles of U, T and q were studied 
for eight time periods. Four levels of data were desired, but signals 
were not always available for all levels due to various component 
failures. Table I summarizes the periods examined and what data were 
available for each period. 
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TABLE I. SUMMARY OF DATA PERIODS 





DATE 


TIME 


NUMBER 

OF 

LEVELS 


NUMBER 

OF 

PROFILES 

(T,U,q) 




1974 








27 


March 


1841-2117 


3 


12 


28 


March 


1902-2059 


3 


6 


13 


August 


1700-1742 


4 


9 


18 


September 


1643 


3 


3 


19 


September 


1735-2243 


3 


24 


20 


September 


0014-0552 


2,3 


15 


21 


November 


1800-2130 


4 


24 




1975 








27 


March 


1502-2106 


2,3 


27 
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Mean wind measurements over ten minute intervals were recorded in 
a log book. The recordings were the number of cup revolutions for the 
interval as recorded by electronic counters. Temperature readings were 
printed on a paper tape from a monitoring digital readout at approxi- 
mately one second intervals. Relative humidity measurements were re- 
corded the same as the temperatures. All information was later coded 
and punched on computer cards for further analysis. 

The data were then screened for gross errors or inconsistencies 

due to instrumentation malfunctions. The criterion at this point for 

retaining or discarding profiles or data periods depended on whether 
2 

values could be computed from temperature spectra for the time 
periods involved. 

The computer analysis of both phases consisted of programming the 
governing equations using either actual data or data picked off the 
profiles. The NPS Fortran program PLOT? was used to plot the predicted 
curves for phases one and two as well as the results from the data. 

1 . Profile Editing 

The data from U, T, and q measurements were plotted on 4-cycle 
semi-logarithmic paper. Since U, T, and q are parameters which vary 
logarithmically with height, a best fit straight line was drawn to the 
data points. In general, the procedures were subjective. In many in- 
stances several different slopes could be obtained from just one graph, 
so a useful criterion was not to give a single data point too much 
influence in determining the straight line. Consequently, the line 
drawn represented a most probable position between data points as 
illustrated in Figure 10. 
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(b) 



(c) 



Figure 10. 






Typical profiles for a) specific humidity - 27 Mar 75 (1539) ; 
b) temperature - 19 Sep 75 (1736) ; and c) winds - 21 Nov 74 
(1842) . 
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Difficulty arose whenever anomalous points appeared in a graph 
as in Figure 11 where the temperature at the third level was obviously 
inconsistent with the second and fourth levels. In a case such as this, 
the point was merely ignored. 

The time period of 21 November 1974 (1817-2130) was particularly 
difficult to analyze due to the apparent discrepancy of the temperature 
and humidity sensors at the fourth level from the readings of the other 
levels. This discrepancy is evident in Figure 12. These sensors were 
found in later tests to be out of calibration, and the decision was 
made to draw the straight line to fit the three lowest levels. 

In some, but not in all, cases the sea surface temperature was 
helpful in determining the slope. However, if the sea surface tempera- 
ture was significantly different from that expected on the basis of the 
temperatures in an atmospheric profile, it could not be used since the 
profile temperatures differed in most cases by only hundredths of a 
degree. 

2 

2. Analysis Procedures to Obtain 

2 

Variance spectra were interpreted to estimate . The spectra 
were obtained using an analog spectral analyzer G Procedures for convert- 
ing temperature spectral values, obtained with the analog spectral 

analyzer, to engineering units and then for obtaining the turbulence 
2 

parameter C T from the temperature spectral estimates are described 
in this sectiono 

a. Spectral Plot Scaling 

A necessary procedure was to scale the spectral plots in 
order to relate RMS input voltages to power spectral densities; (vari- 
ance per unit frequency) . 



40 




Figure 11. Profile of 27 Mar 74 (1929) showing anomalous temperature 
reading at third level. 
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Figure 12. Profiles of 21 Nov 74 (1817). Solid black line represents 
the profile used. 
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To obtain power spectral density levels, corresponding to 
RMS voltage inputs, calibrated scale charts had to be constructed. The 
charts were constructed as follows. Amplitude scaling was accomplished 
by using an externally generated signal whose RMS value was determined 
by a RMS meter. Using a 0-dB ( . Iv) input gain and a spectral gain of 
O-dB (x 1), a signal with amplitude equal to 0.1 volts RMS was produced 
on a screen display as a spike near the selected frequency. This spike 
was then plotted with a X-Y plotter. Successively, the input gain was 
stepped down to attenuate the amplitude of the input signal by 10 dB 
increments, and plots of the height of each resulting spectral spike or 
amplitude was then added to the X— Y plot, at different frequencies. 

These plots represented a graduation of RMS input from a minimum of 
0.001 volts RMS to a maximum of 0.1 volts RMS. An example of such a 
calibration plot on the overlay appears in Figure 13. These procedures 
were performed regularly during analyses to insure continual calibration 
of the spectral analyzer and X-Y plotter. 

For purposes of the chart format the RMS voltages were con- 
verted to LOG 10 units and a graduated scale was constructed so that 
the logarithm of volts RMS could be interpolated from spectral plots. 

The amplitude scale was adjusted for each spectrum as a function of 

both input gain and spectral gain. These values were then converted to 

2 

PSD levels for use in calculating C T values. 

The equation (based on Federal Scientific specifications) 

to convert volts RMS to PSD is 

2 0.79 (cal. level V^) 2 

PSD level (V /Hz) = Filter Bandwidth (2?) 
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Figure 13. Calibration plot 



where the filter bandwidth (3 “ 0.4) is a function of the analysis 
range selected. 

b. Obtaining Turbulence Parameters from Scaled Spectra 

2 

The turbulence parameter was obtained from the tempera- 

ture variance spectra on the basis of the universal formula, Equation 
(2), for the inertial subrange in wave number space S (k) • This expres- 
sion predicts a -5/3 slope for the spectra when plotted in log-log 
formats 0 

Figures 14 and 15 are typical spectra considered in the 
analyses. Temperature spectra often exhibited slopes slightly different 
than the expected -5/3. This feature of temperature spectra have been 
observed by others. The existence of "cold spikes" in temperature 
traces has been given as a reason for this, e.g., Friehe (1976). A 
universally accepted reason is still not available, however. 

Assuming -5/3 slopes for the variance spectra, S (f ) , obtained 
from the analyzer, the intercept of the -5/3 slope with the 1 Hz fre- 
quency line was the spectral density denoted (PSD) value used in comput- 

2 

ing the parameter of interest, . 

The measured PSD value obtained from temperature spectra 
was converted to a spectral density in engineering units by the equation 

S T (f) = C H * PS ° 

= (°C/volt) 2 * Volt 2 /Hz = °C 2 /Hz (28) 

where C is the calibration factor for the temperature system. Determina- 
H 

tion of the value for C has been described by Lund (1975). 

H 

Since the temperature fluctuations were measured at a fixed 
point in the flow, the resultant spectral values are defined at 
"temporal" frequencies, f, and are denoted as S^ff) in equation (28). 
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Figure 14. Typical temperature spectra. ' 2(3) 
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Figure 15 . Anomalous temperature spectra 



To obtain , temporal (f) and space (k) scales had to be 
related in order to use Equation (5) . This was accomplished by using 
Taylor's "frozen turbulence" hypothesis discussed in section II, B. The 
following equation relates temporal to wavenumber spectral values, 



fS„ (f) = kS (k) = c r 2 k" 2/3 
T T 2 T 



(29) 



2 

where S r ^(f) is the spectral density value with units of °C /Hz. This 

2 

leads to the following relationship between C T and S^f), which is 
measured. 



C 



T 



2 




fs T (f) 



(30) 



where, empirically, c^ = 0.25 and k = 2TTf/U. 

From Equation (30), with measured values of f, S^ff), and 

2 

tJ, C can be determined at each level of interest. These were the 
T 

2 

procedures used to determine the variation with height of C T values. 
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IV. RESULTS 



A comparison of the observed results and the Wyngaard et al predic- 
tion curve based on Equation (11) appears in Figure 16 where individual 
data points appear as dots, and averages over Ri intervals of 0.25 
appear as dots within a larger circle. The error bars are standard 
deviations from the mean within each interval, while the number at the 
top of the error bars is the number of observations defining the mean 
value. 

For both the stable (+ Ri) and unstable (- Ri) stratification cases, 
there appears to be little agreement and no noticeable trends. The 
negative Ri numbers were within acceptable ranges from 0.0 to -2.6. The 
positive values, however, ranged from 0.0 to as high as 35.7. Since 
turbulence above a critical Ri number is virtually non-existent due to 
the extreme stability, it was acceptable to assign lower positive Ri 
values to Ri numbers above 1.0. A critical Ri value is generally con- 
sidered to be 0.21. 

The scatter in the observed results in Figure 16 can be attributed 

2 

to scatter in both the measured values as well as 38/3Z values. 

Deviation of temperature spectra from -5/3 slope caused uncertainty in 
2 

estimates. The 90/3Z values are important because the temperature 
gradients are often very small and the quantity is squared for the nor- 
malization. 

In Figure 17, the results computed by Friehe's bulk aerodynamic 
expression. Equation (25) , are compared to the Wyngaard et al prediction 
curve. There seems to be some agreement between the observed results 
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Figure 16. Results of the spectrally derived temperature- structure parameter 
the Richardson number. 




Figure 17. Results of the spectrally derived temperature- structure 
parameter vs. X. 
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and the predicted curve close to near-neutral stability conditions , but 
with the large amount of scatter present , this agreement may be coinci- 
dental. 

The values obtained for Figure 16 and Figure 17 are summarized in 
Tables II and III. 
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Results 



1974 



TABLE II. C 2 and Ri 
T 



Period 


Level 


d0/dz 


°T -3 
x 10 


Ri 




(M) 


°C/cm 




3-27-1841 


4.19 


0.049 


1.67 


0.362 




6.60 


0.031 


3.42 


0.500 




7.64 


0.027 


3.37 


0.500 




13.89 


0.015 


- 


1.000 


3-27-1902 


4.19 


0.047 


2.07 


0.249 




6.60 


0.030 


2.22 


0.392 




7.64 


0.026 


4.34 


0.454 




13.89 


0.014 


- 


0.750 


3-27-1928 


4.19 


0.047 


5.67 


0.230 




6.60 


0.030 


4.98 


0.362 




7.64 


0.026 


6.20 


0.419 




13.89 


0.014 


1.48 


0.750 


3-27-2117 


4.19 


0.019 


3.62 


0.073 




6.60 


0.012 


3.22 


0.115 




7.64 


0.012 


- 


0.133 




13.89 


0.006 


1.75 


0.242 


3-28-1902 


4.19 


0.046 


17.40 


1.000 




6.60 


0.029 


- 


1.500 




7.64 


0.025 


25.90 


1.500 




13.89 


0.014 


1.80 


1.500 


3-28-2059 


4.19 


0.054 


59.60 


0.388 




6.60 


0.034 


- 


0.500 




7.64 


0.030 


30.50 


0.500 




13.89 


0.016 


23.60 


1.000 


8-13-1700 


4.19 


-0.036 


4.55 


-0.538 




6.60 


-0.023 


0.895 


-0.847 




7.64 


-0.020 


1.93 


-0.981 




13.89 


-0.011 


0.281 


-1.783 


8-13-1721 


4.19 


-0.036 


3.95 


-0.665 




6.60 


-0.023 


0.902 


-1.047 




7.64 


-0.020 


2.44 


-1.213 




13.89 


-0.011 


0.338 


-2.204 
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1974 



Period 


Level 


d0/dz 


V -3 


Ri 




(M) 


°C/cm 


x 10 




8 - 13-1742 


4.19 


- 0.025 


4.47 


- 0.451 




6.60 


- 0.016 


0.964 


- 0.711 




7.64 


- 0.013 


2.30 


- 0.823 




13.89 


- 0.007 


0.264 


- 1.496 


9 - 18-1643 


4.19 


0.031 


0.304 


0.308 




6.60 


0.020 


0.161 


0.486 




7.64 


0.017 


- 


0.500 




13.89 


0.009 


0.140 


1.000 


9 - 19-1735 


4.19 


- 0.005 


0.495 


- 0.082 




6.60 


- 0.003 


2.54 


- 0.128 




7.64 


- 0.003 


- 


- 0.149 




13.89 


-0.001 


0.686 


- 0.270 


9 - 19-1846 


4.19 


- 0.005 


5.33 


-0.100 




6.60 


- 0.003 


22.04 


- 0.157 




7.64 


- 0.003 


- 


- 0.181 




13.89 


-0.001 


1.36 


- 0.330 


9 - 19-1900 


4.19 


-0.010 


0.616 


- 0.607 




6.60 


- 0.006 


1.86 


- 0.957 




7.64 


-0 o 005 


- 


- 1.107 




13.89 


- 0.003 


1.12 


- 2.013 


9 - 19-1929 


4.19 


- 0.005 


0.77 


- 0.099 




6.60 


- 0.003 


2.51 


- 0.156 




7.64 


- 0.003 


- 


- 0.180 




13.89 


-0.001 


1.41 


- 0.328 


9 - 19-2011 


4.19 


-0.010 


0.53 


0.051 




6.60 


- 0.006 


1.74 


0.080 




7.64 


- 0.005 


- 


0.093 




13.89 


- 0.003 


0.68 


0.169 


9 - 19-2115 


4.19 


- 0.172 


0.37 


- 3.000 




6.60 


- 0.109 


1.48 


- 3.000 




7.64 


- 0.094 


- 


- 3.000 




13.89 


- 0.052 


0.95 


- 3.000 
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1974 



Period 


Level 


d 0 /dz 


C T -3 


Ri 




( M ) 


° C/cm 


x 10 




9 - 19-2158 


4.19 


- 0.010 


0.48 


- 0.202 




6 . 60 


- 0.006 


2.41 


- 0.319 




7.64 


- 0.005 


- 


- 0.369 




13.89 


- 0.003 


1.35 


- 0.671 


9 - 19-2243 


4.19 


- 0.178 


0.34 


- 3.000 




6.60 


- 0.113 


1.67 


- 3.000 




7.64 


- 0.098 


- 


- 3.000 




13.89 


- 0.054 


1.33 


- 3.000 


9 - 20-0014 


4.19 


- 0.011 


0.69 


- 0.461 




6.60 


- 0.007 


4.55 


- 0.727 




7.64 


- 0.006 


- 


- 0.841 




13.89 


- 0.003 


1.23 


- 1.530 


9 - 20-0056 


4.19 


- 0.013 


0.59 


- 0.408 




6 • 60 


- 0.008 


- 


- 0.643 




7.64 


- 0.007 


- 


- 0.744 




13.89 


- 0.004 


2.32 


- 1.353 


9 - 20-0140 


4.19 


- 0.013 


0.81 


- 0.546 




6.60 


- 0.008 


3.06 


- 0.860 




7.64 


- 0.007 


- 


- 0.995 




13.89 


- 0.004 


1.10 


- 1.809 


9 - 20-0530 


4.19 


- 0.003 


0.80 


- 0.085 




6.60 


- 0.002 


3.79 


- 0.134 




7.64 


- 0.002 


- 


- 0.155 




13.89 


- 0.001 


2.05 


- 0.282 


9 - 20-0552 


4.19 


- 0.006 


0.96 


- 0.952 




6.60 


- 0.004 


- 


- 1.500 




7.64 


- 0.004 


- 


- 1.736 




13.89 


- 0.002 


1.22 


- 3.000 


11 - 21-1800 


4.19 


- 1.626 


3.93 


- 3.000 




6.60 


- 1.032 


19.00 


- 3.000 




7.64 


- 0.892 


7.33 


- 3.000 




13.89 


- 0.490 


7.18 


- 3.000 
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1974 



1975 



Period 


Level 


d0/dz 




Ri 




(M) 


°C/cm 


x T 10 3 




11 - 21-1821 


4.19 


0.008 


3.56 


0.027 




6.60 


0.005 


4.01 


0.043 




7.64 


0.004 


7.14 


0.050 




13.89 


0.002 


2.76 


0.091 


11 - 21-1842 


4.19 


0.006 


3.67 


0.033 




6.60 


0.004 


4.66 


0.051 




7 . 64 


0.004 


7.66 


0.060 




13.89 


0.002 


11.30 


0.108 


11 - 21-1903 


4.19 


0.003 


5.94 


0.009 




6 . 60 


0.002 


6.99 


0.014 




7.64 


0.002 


10.90 


0.017 




13.89 


0.001 


1.53 


0.030 


11 - 21-2027 


4.19 


0.002 


8.14 


- 0.001 




6.60 


0.001 


3.04 


- 0.001 




7.64 


0.001 


13.10 


- 0.001 




13.89 


0.000 


1.10 


- 0.002 


11 - 21-2048 


4.19 


0.002 


8.81 


0.006 




6.60 


0.001 


10.20 


0 o 010 




7.64 


0.001 


12.90 


0.012 




13.89 


0.000 


1.68 


0.021 


11 - 21-2109 


4.19 


- 0.000 


10.20 


- 0.025 




6.60 


-OoOOO 


15.50 


- 0.039 




7. 64 


- 0.000 


14.90 


- 0.046 




13.89 


- 0.000 


1.34 


- 0.083 


11 - 21-2130 


4.19 


0.0 


11.60 


- 0.037 




6.60 


0.0 


17.10 


- 0.059 




7.64 


0.0 


21.40 


- 0.068 




13.89 


0.0 


2.12 


- 0.123 


3 - 27-1503 


4.19 


0.006 


11.10 


0.004 




6.60 


0.004 


13.18 


0.007 




7.64 


0.004 


12.45 


0.008 




13.89 


0 o 002 


- 


0.015 
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1975 



Period 


Level 


dQ/dz 


S 3 


Ri 




(M) 


°C/cm 


x 10 




3 - 27-1519 


4 o 19 


0.006 


10.15 


0 o 004 




6.60 


0.004 


12.66 


0.006 




7.64 


0.004 


5.98 


0.007 




13.89 


0.002 


- 


0.012 


3 - 27-1539 


4.19 


-0.002 


10.55 


- 0.059 




6.60 


-0.001 


12.57 


- 0.092 




7.64 


-0.001 


9.36 


- 0.107 




13.89 


- 0.000 


- 


- 0.195 


3 - 27-1600 


4.19 


- 0.005 


11.25 


- 0.032 




6.60 


- 0.003 


11.63 


- 0.051 




7.64 


- 0.003 


9.63 


- 0.059 




13.89 


-0.001 


- 


- 0.107 


3 - 27-1644 


4.19 


-0.002 


11.82 


- 0.014 




6.60 


-0.001 


16.01 


- 0.023 




7.64 


-0.001 


10.92 


- 0.026 




13.89 


- 0.000 


- 


- 0.048 


3 - 27-1656 


4.19 


0.0 


11.88 


- 0.013 




6.60 


0.0 


12.33 


-0.020 




7.64 


OoO 


7.56 


- 0.024 




13.89 


0.0 


- 


- 0.043 


3 - 27-2020 


4.19 


0 o 041 


- 


0.016 




6.60 


0.026 


- 


0.025 




7.64 


0.022 


10.02 


0.029 




13.89 


0.012 


- 


0.053 


3 - 27-2045 


4.19 


0 o 044 


5.43 


0.025 




6.60 


0.028 


- 


0.039 




7.64 


0.024 


0.72 


0.045 




13.89 


0.013 


- 


0.082 


3 - 27-2106 


4.19 


0.015 


5.73 


- 0.046 




6.60 


0.009 


- 


- 0.072 




7.64 


0.008 


7 0 41 


- 0.083 




13.89 


0 o 004 


- 


- 0.151 
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1974 



TABLE III. C,^ and X Results 



Period 


Level 


(T - T ) 


C 


X 




(M) 


a °c s 


x T lcf 3 




3-27-1841 


4.19 


2.76 


1.67 


_ 




6.60 


3.06 


3.42 


0.3121 




7.64 


3.18 


3.37 


0.3688 




13.89 


3.22 


- 


- 


3-27-1902 


4.19 


2.89 


2.07 


_ 




6.60 


3.23 


2.22 


0.2796 




7.64 


3.32 


4.34 


0.3266 




13.89 


3.32 


- 


- 


3-27-1928 


4 o 19 


2.97 


5.67 


_ 




6 o 60 


3.30 


4 o 98 


0.1966 




7.64 


3.54 


6.20 


0.2381 




13.89 


3.45 


1.48 


0.4268 


3-27-2117 


4.19 


3.25 


3.62 


_ 




6 . 60 


3.57 


3.22 


0.2418 




7.64 


3.59 


- 


- 




13.89 


3.63 


1 0 75 


0.4863 


3-28-1902 


4.19 


0.86 


17.40 


-0.3101 




6.60 


0.99 


- 


- 




7.64 


0.98 


25.90 


0.3046 




13.89 


1.32 


1.80 


0.7809 


3-28-2059 


4.19 


1.99 


59.60 


-0.0386 




6.60 


2.36 


- 


- 




7.64 


2.37 


30.50 


0.4924 




13.89 


2.60 


23.60 


1.0340 


8-13-1700 


4.19 


-0.59 


4.55 


-0.1676 




6.60 


0.69 


0.895 


0.0849 




7 o 64 


0.49 


1.93 


0.0375 




13.89 


0.38 


0.281 


0.0013 


8-13-1721 


4.19 


-0.70 


3.95 


-0.1893 




6 0 60 


0.64 


0.902 


0.0691 




7.64 


0.43 


2.44 


0.0223 




13.89 


0.39 


0.338 


0.0051 
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X 



1974 



Period Level 

(M) 




S 2 

x 10 



-3 



8 - 13-1742 


4 o 19 


i 

o 

• 

CO 

-0 


4.47 


- 0.2270 




6*60 


0.34 


0.964 


- 0.0054 




7.64 


0.18 


2.30 


- 0.0435 




13.89 


0.16 


0.264 


- 0.1194 


9 - 18-1643 


4.19 


0.24 


0.304 


0.0229 




6.60 


0 o 19 


0 o 161 


- 0.0523 




7.64 


0.22 


- 


- 




13.89 


0.41 


0.140 


0 O 0443 


9 - 19-1734 


4 0 19 


- 0.47 


0 o 495 


- 0.0821 




6.60 


- 0.53 


2.54 


- 0.1554 




7.64 


- 0.52 


- 


- 




13.89 


- 0.58 


0.686 


- 0.3492 


9 - 19-1846 


4.19 


- 0.66 


5.33 


- 0.1453 




6 o 60 


- 0.73 


22 o 04 


- 0.2409 




7.64 


- 0.70 


- 


- 




13.89 


- 0.81 


1.36 


- 0.6136 


9 - 19-1900 


4.19 


- 0.60 


0.616 


- 0.1542 




6.60 


- 0.65 


1.86 


- 0.2459 




7.64 


- 0.62 


- 


- 




13.89 


- 0.77 


1.12 


- 0.5920 


9 - 19-1929 


4.19 


- 0.70 


0.77 


- 0.1395 




6.60 


- 0.74 


2.51 


- 0.2209 




7 0 64 


- 0.66 


- 


- 




13.89 


- 0.81 


1.41 


- 0.4974 


9 - 19-2011 


4.19 


- 0.62 


0.53 


- 0.1174 




6.60 


- 0.67 


1.74 


- 0.1848 




7.64 


- 0 o 61 


- 


- 




13.89 


- 0.78 


0.68 


- 0.4478 


9 - 19-2115 


4.19 


- 0.40 


0.37 


- 0.1036 




6.60 


- 0.42 


1.48 


- 0.1515 




7.64 


- 0.37 


- 


- 




13.89 


- 0.55 


0.95 


- 0.3878 
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Period 


Level 

( M ) 


(T 


9 - 19-2158 


4.19 


- 0 , 




6.60 


-0 




7.64 


-0 




13.89 


-0 


9 - 19-2243 


4.19 


-0 




6.60 


-0 




7.64 


-0 




13.89 


-0 


9 - 20-0014 


4.19 


-1 




6.60 


-1 




7.64 


-1 




13.89 


-1 


9 - 20-0056 


4.19 


-1 




6.60 


-1 




7.64 


-1 




13.89 


-1 


9 - 20-0140 


4.19 


-1 




6.60 


-1 




7.64 


-1 




13.89 


-1 


9 - 20-0530 


4.19 


-1 




6.60 


-1 




7.64 


-1 




13.89 


-1 


9 - 20-0552 


4.19 


-1 




6.60 


-1 




7.64 


-1 




13.89 


-1 


11 - 21-1800 


4.19 


0 




6.60 


0 




7.64 


0 




13.89 


0 



C T 2 -3 
x 10 


X 


0.48 


- 0.0673 


2.41 


- 0.0987 


1.35 


- 0.2496 


0.34 


- 0.1702 


1.67 


- 0.2963 


1.33 


- 0.6489 


0.69 


- 0.3783 


4.55 


- 0.5992 


1.23 


- lo 3172 


0.59 


- 1.4922 


2.32 


- 4.4630 


0.81 


- 1.5780 


3.06 


- 2.1809 


1.10 


- 5.0870 


0.80 


- 0 o 4668 


3.79 


- 0.7335 


2.05 


- 1.5530 


0.96 


- 0.6111 


1.22 


- 2.2413 


3.93 


0.0107 


19.00 


0.0393 


7.33 


0.0552 


7.18 


0.0480 



- T 

°C s 

42 

45 

38 

56 

67 

88 

76 

91 

37 

48 

44 

55 

56 

68 

69 

74 

54 

67 

66 

73 

60 

63 

61 

67 

69 

72 

69 

82 

54 

57 

63 

,58 



60 



X 



1974 



1975 



Period 


Level 


(T - T ) 


2 

c 




( M ) 


d °c s 


T 

x 10 


11 - 21-1821 


4.19 


0.51 


3.56 




6 . 60 


0.53 


4.01 




7.64 


0.57 


7.14 




13.89 


0.53 


2.76 


11 - 21-1842 


4.19 


0.45 


3.67 




6.60 


0.47 


4.66 




7.64 


0.52 


7.66 




13.89 


0.46 


11.30 


11 - 21-1903 


4.19 


0.38 


5.94 




6.60 


0.39 


6.99 




7.64 


0.45 


10.90 




13 0 89 


0.38 


1.53 


11 - 21-2027 


4.19 


0.05 


8.14 




6.60 


0.04 


3.04 




7.64 


0.10 


13.10 




13.89 


0.10 


1.10 


11 - 21-2048 


4.19 


- 0.01 


8.81 




6.60 


- 0.01 


10.20 




7.64 


0.05 


12.90 




13.89 


0.0 


1.68 


11 - 21-2109 


4.19 


- 0.03 


10.20 




6.60 


- 0.05 


15.50 




7.64 


0.0 


14.90 




13.89 


- 0.08 


1.34 


11 - 21-2130 


4.19 


- 0.01 


11.60 




6.60 


- 0.02 


17.10 




7.64 


0.04 


21.40 




13.89 


- 0.06 


2.12 


3 - 27-1502 


4.19 


- 0.78 


11.10 




6.60 


- 0.70 


13.18 




7.64 


- 0.77 


12.45 




13.89 


- 0.77 


- 



0.0136 

0.0373 

0.0591 

0.0452 

0.0052 

0.0242 

0.0474 

0.0208 

- 0 o 0046 
0.0044 
0.0260 
- 0.0125 

- 0.0729 

- 0.0848 

- 0.0760 

- 0.2260 

- 0.1442 

- 0.1637 

- 0.1591 



- 0.1697 

- 0.2066 

- 0.5626 

- 0.1328 

- 0.1572 

- 0.1596 

- 0.4471 

- 0.1804 

- 0.1056 

- 0.1394 



61 



X 



1975 



Period Level (T - T ) C 





(M) 


a *c s 


T - J 
x 10 




3-27-1519 


4.19 


-0.18 


10.15 


-0.1450 




6 • 60 


-0.35 


12.66 


-0.0736 




7 o 64 


-0.38 


5 0 98 


-0.1009 




13.89 


-0.34 


- 


- 


3-27-1539 


4.19 


-0.74 


10.55 


-0.1528 




6.60 


-0.71 


12.57 


-0.0924 




7.64 


-0.72 


9.36 


-0.1164 




13.89 


-0.75 


- 


- 


3-27-1600 


4.19 


-0.73 


11.25 


-0.1400 




6.60 


-0.69 


11.63 


-0.0836 




7.64 


-0.71 


9.63 


-0.1035 




13.89 


-0.76 


- 


- 


3-27-1644 


4.19 


-0.74 


11.82 


-0.1164 




6.60 


-0.62 


16.01 


-0.0615 




7.64 


-0.71 


10 o 92 


-0.0825 




13.89 


-0.74 


- 


- 


3-27-1656 


4.19 


-0.32 


11.88 


-0.0838 




6.60 


-0.43 


12.33 


-0.0438 




7.64 


-0.42 


7.56 


-0.0537 




13.89 


-0.46 


- 


— 


3-27-2020 


4.19 


0.85 


- 


- 




6.60 


0.72 


- 


- 




7.64 


1 o 08 


10.02 


0 o 0061 




13.89 


1.08 


- 


— 


3-27-2045 


4.19 


0.15 


5.43 


-0.3261 




6.60 


0.50 


- 


- 




7.64 


0.60 


0.72 


-0.1205 




13.89 


0.63 


- 


— 


3-27-2106 


4.19 


0.34 


5.73 


-0.2351 




6.60 


0.44 


- 


- 




7.64 


0.40 


7.41 


-0.1566 




13.89 


-0.24 


- 


- 
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V. CONCLUSIONS 



2 

Correlation of spectrally derived C T results with the stability 

parameter z/L or with the Richardson number values was not good for 

these data in a marine environment. The temperature gradient in a 

2 

marine environment is too small and gives rise to very large C^ values 
which contribute to the observed scatter. As pointed out by Friehe 
(1976) , a complication in the marine boundary layer is the sometimes 
anomalous behavior of the temperature field where inertial subranges 
are not always obtained. Salt spray is also a factor to content with 
in a marine environment and can be the cause of many problems with 
instruments exposed to it. 

For Navy operational purposes, then, existing empirically derived 
expressions defined by z/L or the Richardson number may require modi- 
fication before use in describing the atmosphere over an ocean environ- 

2 

ment for obtaining C T indirectly. However, for reasons explained m 

the preceding paragraph, more measurements using a different approach 

2 2 
in computing may lead to better correlation between C T and an 

easily defined stability parameter. It should also be noted that the 

correct empirical relationship, when finally derived, may be different 

over the marine environment than the curve predicted by Wyngaard et al 

over land. 
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